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_ - A class of random graphs is introduced and studied. The graphs are constructed in an algorithmic way from 

' five motifs which were found in [Mile R., Shen-Orr S., Itzkovitz S., Kashtan N., Chklovskii D., Alon U., Science, 

2002, 298, 824-827]. The construction scheme resembles that used in [Hinczewski M., A. Nihat Berker, Phys. 
^ I Rev E, 2006, 73, 0661 26], according to which the short-range bonds are non-random, whereas the long-range 

^ , bonds appear independently with the same probability A number of structural properties of the graphs have 

^ been described, among which there are degree distributions, clustering, amenability, small-world property For 

one of the motifs, the critical point of the Ising model defined on the corresponding graph has been studied. 
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1. Introduction and setup 



■3 

. A vast variety of large systems occurring in nature and society have a very complicated topolog- 

ical structure. These are the Internet, the World Wide Web, citation, neural, and social networks, 
etc. In view of the complex topology and unknown organizing principles, the networks are often 
modeled as random graphs. A random graph with a given node set y is a graph in which for a given 
^ . pair i,j S V, the bond («, j) appears at random. The study of random graphs has been originated 

' by P. Erdos and A. Renyi, who were the first to introduce such graphs in In the Erdos-Renyi 

0^ ■ random graph model, denoted by Gn,p , the number of nodes is n, and the bonds between distinct 

' nodes appear independentljQ with the same probability p. An important characteristic of a graph 

is the node degree, which is the number of bonds attached thereat. In Gn,p , it is a random variable. 



' and all such variables are independent and have the same probability distribution. Namely, the 



probability that the degree of a given node is k is given by the Bernoulli law 

p„(fc) = (""^ fc = o,i,...,n-i. (1) 



For random graphs, the most interesting questions refer to their asymptotic properties in the limit 
' n — +O0. To get nontrivial answers to such questions one allows the parameters to depend on n. 

For pn — c/n, the limit of ([TJ is the Poisson law 

P(fc) = c'^'e^Vfe! , fee No. (2) 

However, for irregular complex networks, the random graph Gn,p is not a good model since in 
the most of such networks the degree distribution is essentially non-Poissonian. Many real world 
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networks, e.g. the WWW, are characterized by power-law node degree distributions, which have 
the form P{k) = Ck~^, 7 > 1, typical for the so-called scale-free graphs. 

Another important parameter characterizing a random graph is the clustering coefficient, which 
is the probability that two nodes are neighbors given they have a common neighbor. Clearly, in 
Gn,p this probability is p, and is the same independently of whether or not the nodes have a com- 
mon neighbor. Real world networks usually manifest strong clustering, which once more indicates 
that Erdos-Renyi type random graphs are not appropriate as their models. In Q, D.J. Watts and 
H. Strogatz proposed another type of random graphs, in which the mentioned disadvantage is 
overcome. It should be noted, however, that such graphs do not have power law node degree distri- 
butions. The next step beyond the Erdos-Renyi model was done by A.-L. Barabasi and R. Albert 
in ||5j. In their model, the preferential attachment principle has been employed, typical for many 
real networks (the Internet, citation and social networks). According to this principle, the more 
connected a node is, the more likely it receives a new bond. The construction of the Barabasi- 
Albert model starts from an initial graph with mg ^ 2 nodes (neither can be isolated). At each 
step, one adds a new node and connects it to the existing nodes. The probability pi that the new 
node is connected to node i is proportional to the degree rii of that node, that is, pi — ni/(^j nj). 

Many of the complex networks occurring in nature contain characteristic patterns, recurring 
much more frequently than the other ones. They are called network motifs, see (6-ll|. Different 
networks may have different motifs, and motifs in turn can characterize the networks. For instance, 
in biological regulation networks it has been experimentally demonstrated that each of the mo- 
tifs can perform a key information processing function, see In [l^, the authors introduced a 
random graph model based on some geometric principles (constraints). Then, they compared the 
appearance of eight elementary three- and four-node patterns in their model with the same charac- 
teristics of the Erdos-Renyi random graph. It turned out that five of these patterns are motifs for 
their model, but not for the Erdos-Renyi random graph, see figure [T] below and table 1 in fld\. For 



Figure 1. Three and four node motifs Mi, M2, M3, M4, M5 found in [T3|. 

another random graph model, the appearance of the same eight patterns was also studied in 
Note that among these patterns, only Mi and M5 correspond to complete graphs (each node is a 
neighbor to every other node). 

One of the ways to get information about infinite graphs, also random ones, is to study the 
properties of certain models of statistical physics defined thereon. The most po pula r ones are the 
Ising and Potts models, as well as the models of bond and site percolation, see |l2| . On the other 
hand, in statistical physics certain graphs are employed to mimic a crystal lattice, for which the 
critical-point behavior of the Ising model can be described in an explicit and rigorous way. These are 
the so-called hierarchical lattices introduced in jlsl . Such lattices are defined in a rather algorithmic 
way by means of the basic pattern, e.g., by a "diamond", which is the pattern M3 depicted in figure[T] 
A mathematical description of the Gibbs states of the Ising model on such graphs was done by 
RM. Bleher and E. Zalys in M. Hinczewski and A. Nihat Berker ^ studied the critical- 

point properties of the Ising model on the diamond hierarchical lattice 'decorated' by additional 
bonds, which appear at random. In the present paper, we follow the way suggested in (l6| and 
introduce hierarchical graphs constructed by means of the motifs shown in figure [T] decorated 
by additional bonds which appear at random and repeat, in a way, the corresponding motif. We 
analyse some of their characteristics, such as the average degree, the node degree distribution, 
amenability, the small-world property, as well as the critical-point properties of the Ising model. 
This study of ours is a continuation of jl3|, where the graph based on Mi was introduced by one 
of us. Note that our hierarchical graphs also found applications in cryptography, see |l8| . 
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2. The graphs: construction and structural properties 

2.1. The construction: informal description 

As is typical for hierarchical graphs, e.g., for hierarchical lattices in fl3l[l6t. the construction is 
carried out in an algorithmic way: at fc-th level, fc g N, one produces a subgraph, say , which is 
then used as a construction element for producing A^+i . The procedure is the same at each level, 
and the starting element is obtained from the corresponding motif. Let us illustrate this for the 
simplest case based on the motif A/i . First we label the nodes of Mi by a, b, and c, as it is shown 
in figure and obtain Ai — the triangle. Then we take three such graphs and label them by the 
same labels. As a result, each of the triangles has nodes of two kinds: one node the label of which 
coincides with the triangle label (e.g., node a in triangle a), and two nodes with non-coinciding 
labels. Thereafter, the triangles are glued up according to the following rule: node c of triangle a 
is glued up with node a of triangle c, ect. The nodes with the coinciding labels remain untouched. 
These are the so-called external nodes of A2 . The remaining nodes are called internal. The bonds 
of the initial triangles turn into the bonds of A2 . We call them basic bonds; they are depicted by 
solid lines. At the next stage of the first step, we add the bonds connecting the external nodes in 
the same way as it is in the motif Mi . Such bonds are depicted by dotted lines and are called 
decorations. As a result, we obtain the graph A2 , which has six basic bonds and three decorations, 
three external and three internal nodes. Then we repeat the same procedure — take three copies 
of A2 , label them by a, 6, and c, and divide their external nine nodes into two groups: three nodes 
with coinciding labels and six nodes with non-coinciding labels. Then the graphs A2 are glued up 
as described above. Thereafter, three decorating bonds are drawn to connect the external nodes. 
This procedure is repeated ad infinitum. Similar constructions for motifs A/2 and M3 are presented 
in figure |3] and figure H) respectively. Here we have omitted the decorating bonds not to overload 
the pictures. The picture for M4 is obtained from that for M^ by adding the diagonals. The picture 
for M5 is just the three-dimensional version of the picture for Mi , where the basic pattern is a 
tetrahedron. 

2.2. Definitions 

In order to fix the terminology and to make the construction of our graphs mathematically 
immaculate, we begin by introducing a number of relevant mathematical notions. A (simple) graph 
G is a pair of sets (V, E), where V is the set of nodes, whereas E is a subset of the Cartesian product 

V X V. We suppose that E is symmetric and irreflexive, i.e., (j, i) e E whenever {i,j) e E, and 
(i, i) ^ E for every i,j G V. We say that i and j are connected by a bond if {i,j) £ E. In this case, 
we write « ~ j and say that i and j are adjacent or that they are neighbors. Hence, the elements of 
E themselves can be called bonds. The graph is said to be complete if each two nodes are adjacent. 
For a given i, by n{i) we denote the degree of i — the number of its neighbors. If V is finite, then E 
is also finite, and the graph is said to be finite. Otherwise, the graph is infinite. An infinite graph 
is called locally finite, if n(i) is finite for every node. All the infinite graphs we study are countable, 
which means that both sets V and E are infinite and countable. Given G = (V, E) and G' = (V, E'), 
let : V — > V be such that 4>{i) <t>{i) whenever i j. Such a map (p is called a morphism. A 
bijective morphism is called an isomorphism. If 4> is an isomorphism, then its inverse (f>^^ is also 
an isomorphism, and then the graphs G and G' are said to be mutually isomorphic. Such graphs 
have identical structures and thus can be identified. In this case, we also say that G' is a copy of 
G. One observes that this refers to both finite and infinite graphs. An isomorphism : V — > V, i.e. 
which maps the graph onto itself, is called an automorphism. The graph G' = (V, E') such that 

V C V and E' C E is said to be a subgraph of G = (V, E). In this case, we write G' C G. Suppose 
that a subgraph G' C G has a copy, say G, that is, there exists an isomorphism cf) : Q ^ G . Then 
(p, considered as a map ^ : G G, is called an embedding of G into G, whereas G' is called the 
image of G under this embedding. Figure [T] presents the so-called unlabeled graphs, which we call 
patterns. After labeling, i.e., attaching a label to each of the nodes, such a pattern turns into a 
graph. Another labeling may or may not yield the same graph. This depends on whether or not 
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there exists the corresponding automorphism. For instance, any labeHng of the triangle Mi yields 
the same graph since in any case each of the nodes has the same neighbors. For the pattern M2 , 
the left-hand graph in figure[3]with the interchanged labels a and b is the same. However, the graph 
with the interchanges c and d is not the same anymore. Of course, this new graph is isomorphic 
to the initial one. This is because there is only one nontrivial automorphism of M2 : the one which 
interchanges a and b, and preserves c and d. The triangle has six automorphisms. 

Let us now turn to random graphs. To introduce such a graph we need an underlying graph 
G = (V, E) and a family £ of subsets of E. If G is finite, as £ one can take the set of all subsets 
of E. In the sequel, we deal with such graphs only. Thus, for E' £ £, we say that E' has been 
picked at random with probability /'(E'). In the Erdos-Renyi model Gn,p , the underlying graph 
is complete with the node set V = {1, . . . ,n}. In this model, P(E') = p'^ where p € [0, 1] and 
|E'| stands for the number of elements in E'. In other words, the elements of E are being picked 
independently, each with the same probability p. In a bit complicated model, the bonds are picked 
independently but with probability which depends on the bond. In this case, as well as in the case 
of the Erdos-Renyi model, we deal with a random graph with independent bonds. For such graphs, 

p(E')- n^'(^)' (3) 

eSE' 

where p{e) is the probability of picking bond e. The set of graphs G' = (V, E') with E' E £ is called 
the graph ensemble — each G' is picked at random from this ensemble. A random graph model is 
the pair consisting of the graph ensemble and of the function P : £ [0,1]. If the function P is 
as in ([3]), the graph is said to be a random graph with independent bonds. Suppose that we have 
two random graph models with independent bonds. Let also G = (V, E) and G = (V, E) be their 
underlying graphs and : V — )■ V be a morphism. Then this map is said to be the morphism of the 
random graphs if for every £ E, the probability (in the first model) that this bond is picked 

is the same as the corresponding probability (in the second model) for the bond {(pii) , <j){j)) . 

2.3. The construction 

As was mentioned above, each of our graphs is constructed in an algorithmic way from the cor- 
responding motif presented in figure [T] Since they are going to be random graphs with independent 
bonds, we have to construct the corresponding underlying graphs and, to define, for a given bond, 
the probability of being picked, c.f. (|3]). In all our models, the bonds will be of two kinds, which 
we call basic bonds and decorations. Basic bonds are going to be non-random, i.e. they are picked 
with probability one. Decorating bonds appear with probability p £ [0, 1], which is a parameter of 
the model. Turn now to the construction of the underlying graphs. Let q be the number of nodes 
in the corresponding motif, that is, 9 = 3 for Mi and q = A for the remaining motifs. At step 
fc = 1, we just label the vertices of the corresponding motif by i = 1, . . . , q and obtain the initial 
graph Ai = (Vi ,Ei). All its bonds are set to be basic. Suppose now that we have q + 1 copies 
of Ai obtained by the isomorphisms (/)2 , j = 0, 1, . . . , q. Thus, in j-th copy the nodes are 9!>2(*), 
i = 1, . . . ,q. The graph A2 is obtained from these copies under the following conditions 

'/'2(*) = <^2(*), i = l,...,q; 4>\{j) ^ i = l,...,q, i^j. (4) 

Thus, the images of Vi under 02 and 02 with i ^ j intersect only at one node where (U} holds. The 
maps (1)2 , j = 0,1, . . . , q embed Ai into A2 . The nodes 02(*); * — 1: ■ • ■ ? 9j are called the external 
nodes of A2 . All other nodes are called internal. Thus, A2 has q external and q{q — l)/2 internal 
nodes. At this stage, we label them by i = 1, . . . , q(q + l)/2 in such a way that the external nodes 
have the same labels as in Ai , that is, 02(0 = *i * — Ij • ■ • 1- -By construction, the bonds obtained 
as images under the map 0° decorations: they are of the form (02(*)i </^2(j)) where i and j 
are adjacent in Ai . From the first condition in (|4]) we see that the decorating bonds connect the 
external nodes of A2 . The remaining bonds of A2 are set to be basic. Now we construct A3 from 
one copy of Ai and q copies of A2 . Let 0° be the map which produces the copy of Ai and 03 , 
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j = 1, . . . ,qhe the maps which produce the copies of A2 . We then impose the conditions 

and obtain A3. Thus, embeds Ai — ^ A3 , and : A2 A3 . « = 1, 2, . . . , q. As above, the nodes 
4>^^{i) are set to be external, and the remaining nodes are internal. The images of V2 under (jj^^ and 
(/)g with i j intersect only at one node where ^ holds. Again we label the nodes of A3 in such 
a way that </>|(j) = i, i — 1,. . . ,q. Now let us establish which bonds of A2 are decorating and 
which are basic. As above, the bonds connecting the external nodes are decorating. The images of 
decorating bonds of A2 are decorating bonds in A3 ; the same is also true for the basic bonds — the 
basic bonds of A3 are exactly the images of the basic bonds of A2 . For fc ^ 4, the construction of 
Afc from Afe_i and Ai is identical to the construction of A3 just described. As above, by Vk and Ek 
we denote the sets of nodes and bonds of A^ , respectively. Thus, for k ^ 2 we have Ek = E'^VJ E'^ , 
where E'^, (respectively, E'^) consists of basic (respectively, decorating) bonds. All Afc , fc G N, are 
considered as subgraphs of an infinite graph Aqo , the structure and properties of which are not 
important for the study presented in this article. 

Note that the construction principle used above essentially differs from that used in the con- 
struction of the hierarchical lattices in flsl - flGj . Namely, in our case to obtain A^ one replaces each 
node of the basic pattern by a copy of the graph Afc_i . In the hierarchical lattices, one replaces 
a bond. As we shall see in the sequel, this leads to essentially different properties of the resulting 
graphs. Below in figure |2] we illustrate the construction described above for the case where the 
basic pattern is the motif Mi . In this case, the bare graph (i.e. the one which occurs for p = 0) 
is the approximating graph for the fractal known as the Sierpihski gaskei@. The elements of E2 
(middle graph) and of £"3 (right-hand graph) are depicted by solid lines. The elements of E2 and 
of i?3 are depicted by dotted lines. We omit some dotted lines to indicate that they appear at 
random and hence may be absent in a given realization of the graph. Note that A3 can be viewed 
as the triangle composed of three copies of A2 . In figure [31 we present the construction of the 
bare graph A3 corresponding to M2 ■ In contrast to the former case, this is not a planar graph. 
In figure U] we construct the bare graph A2 for the motif A/3 . One observes that in that picture 
the node c of the lower left-hand quadrat (i.e. quadrat a) is glued up with node a of the upper 
right-hand quadrat. It is interesting that the corresponding fractal can be obtained by the following 
procedure, resembling the one which yields the Sierpihski gasket. One takes the full quadrat and 
cuts it out into four equal quadrats, but without cutting the external lines. Then, one glues up the 
vertices of the smaller quadrats as depicted and proceeds with cutting out the smaller quadrats. 
The fractal which one obtains from M5 is a three dimensional version of the Sierpihski gasket. One 
takes the full tetrahedron and cuts out its inner one fourth in such a way that the remaining four 
tetrahedra are glued up according to the rule: vertex b of tetrahedron a is glued up with vertex a 
of tetrahedron b, etc. 



c 




Figure 2. Construction of the graph A3 based on Mi . 



^The fractal itself is obtained as the closure of the set U^gpsiVfe in the appropriate topology, see e.g. |19( . 
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Figure 3. Construction of the bare graph A3 based on M2 




Figure 4. Construction of the bare graph A2 based on M3 . 



2.4. Degree distribution 

Now we turn to the description of the structural properties of the graphs constructed above. Let 
mfc , fc e N, be the number of times the basic pattern appears in non-decorated Afc as a subgraph. 
For the graphs based 011 AIi and Mr-, we have the same situation. Here mi = 1 and ruk = qmk-i , 
for A; ^ 2 with the exception in A2 , where the additional pattern appears. So, for Mi and M5 we 
obtain mi = 1 and 

mfe = + k^2. (6) 

Here q = 3 and = 4 for Mi and M5 , respectively. For motif M3 we have mi = 1, m2 = 2g and 

mk = qmk-i, > 3. (7) 

Hence irik = 2 • for k ^ 2. The simplest case if for M4 , where nik = qm,k-i for A; ^ 2. It 
gives m4 = 4J^~^. The last motif is M2 - triangle with additional bond. On each level this bond 
"produces" new 17 patterns. So mi = 1 and for fc > 2 ruk = 5(26 • 4*^~^ — 17). 

Now we analyze the number of times the basic pattern appears in fully decorated Afc , denoted 
by ffik ■ For the graphs based on Mi and M5 wc have 



mfe = 



2q+l 
1^' 



q-1 



(8) 



For M3 it is 14*^ — | and for M4 we obtain 1(4*^ — !)• In all cases, we have mfc increasing as 
C{p)q'^~^, which means that adding decorations does not change the asymptotics of mfe . 
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In a similar way, we obtain 

\Vk\^^^, |i?fc|=rg^-i+rp^^^, fc G N, (9) 

where \Vk\ stands for the number of nodes in , whereas \Ek\ is the expected number of bonds 
in this graph. 

As was mentioned above, the degree distribution is a very important characteristic of the graph. 
In contrast to Erdos-Renyi type graphs, for our graphs the distribution of the random variable n{i) 
depends on the type of i. Thus, the simplest way to describe this distribution is to average n{i) 
over the nodes of a given A^. , that is, to consider 

Let (n/c) be the expected value of rife . Then 

= ^E.imi - ^('^ - 1 - ^ • '-^^ ■ (11) 

However, this result provides only partial information on the node degree distribution. To get more, 
let us analyze the structure of the node sets Vk , k — 1,2, . . . , more in detail. For a given A^ and 
^ = 1, . . . , fc, let V^''^ be the set of nodes i G Vk which are external for some A/ and, at the same 
time, are internal for any A/_|_i . Of course, here we mean those A/'s which are subgraphs for A^ . 
As an example, let us consider the graph A2 based on Mi , see the middle graph in figure [21 The 
nodes a, b, and c constitute V2 '^\ whereas the remaining nodes constitute V2 ^\ 

The elements of V^'^ are exactly the nodes at which the subgraphs A;^ , j — 1,. . . ,q are 

(k) 

glued up to form A^ , whereas the elements of V^: are exactly the external nodes of A^ . Then 

\V^'''^ \ = q and \V^''~'^'^\ = q{q - l)/2. For « < fc - 1, we have \V^^'' \ = q\V^'2i\, which can be solved 
to yield 



\V^'^\ = 7:<l'-Hq-l), l^l,...,k-l, \Vl'^\=q. (12) 



The reason to consider the sets V^^^ is that all the elements of each such V^^^ have the same 
degree distribution, independent of A; for Z ^ fc — 1. In fact, the degrees of i e 14^^^ are non-random 
since these nodes receive no decorating bonds. For such i, n{i) = 'n'^^Hj)^ where n'^°^(j) is the 
degree of the corresponding node in the basic pattern, and the sum is taken over all such patterns 
which are glued up. By the symmetry of Mi , M3 , and il/5 , we have n(i) = 4 for Mi and il/3 , and 
n{i) = 6 for M5 . For M2 , n{i) takes values 3, 4, 5, see figure[3] For i € V^''^ , I = 2,3, . . . ,k — 1, we 
have n{i) = n(i) + v{i), where n{i) is non-random and has to be calculated as just described. The 
summand v{i) is the number of decorating bonds attached to i. To simplify our consideration, let 
us stick to the case of the graph generated by Mi . Then for ^ = 1, . . . , fc — 1 and i e vj:'"'^ , we have 
n{i) ~ 4 and v{i) takes values 1/ = 0, 1, 2, . . . , 4(Z — 1), with probability 

Vmh{v{i) = v)= (^^^^ -P)*^'"'^"". (13) 

For i G y^^^ , takes values 0, 1, ... , 2(fc — 1). Therefore, the maximum node degree which can 
occur in Vfc , fc ^ 2, is 

maxn(i) = 4A; - 4. (14) 

As is usual in the theory of real world networks, which are in fact non-random, the randomness 
manifests itself as the random choice of a node. If we apply this principle here, then (|13p can be 
considered as the conditional probability distribution, conditioned by the event that the node i 
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has been picked from the set V^^\ The probabiHty of the latter event is taken to be proportional 
to the number of its elements, that is, 

Prob 



Prob 



^ 



2q 



If we now take the expectation of n{i) with respect to this distribution, that is, 

1 = 1 u=0 ^ V / 

2(fe-l) 



E (2 + '^)^.^•('^"'Va-p)^^'=-^^-^ 

z>— ^ ^ 



we readily obtain (nfe) = 4 + 2p — 4(l-|-p)/ (3'^"^ + 1) , which is in full agreement with the one given 
in (fTTj) or in table [T] In order to figure out the limit k — >■ +oo of the distribution given by (p!3|l 
and (|15p we calculate its characteristic function, c.f. (|16p . 

^fcW = E E ^-P[i*(4 + -)]TTi^ ■ 'Mp^(l-P)^('-^'-'' (17) 

1=1 u=0 \ ^ J 

2(fe-l) 



E exp[it(2 + .)]^^^ .(2(^;i))p^(i-p)^('-)-^ 

26^'* 1 - 31-'= (e"39+ 1 -p)^ 



l + 31-fe 3 _ (ei*p + 1 - 

96^" 2 

Then the limiting characteristic function is 

2e^'* 

3 — {e}^p + 1 — p) 

which can be continued to a meromorphic function analytic in some complex neighborhood of 
the point t — Q. This means that the limiting node degree distribution has all moments and 
hence cannot be of scale-free typ^. This also agrees with the heuristic rule, see equation (3.13) on 
page 188 in that for scale-free networks 

maxn(i) - |Vfc|i/("^i', a > 1, (19) 

whereas in our case we have (HH) and 1141 = (3*= + 3)/2. In a similar way, one can show that all 
our graphs are not scale-free. Another observation on this item can be made by comparing the 
function (|18p with the characteristic function of the Poisson distribution ^ which has the form 

'^Poisson(0 = exp [c (e'* - 1)] , 

and hence can be continued to a function analytic on the whole complex plane. Therefore, the 
distribution corresponding to psp with p > is intermediate as compared to the Poisson and 
scale-free distributions. For p — Q, the function (fT5)) is also entire. 



^For scale-free graphs, the node degree distribution is P{k) = Ck ^, fc 1, 7 > 1; hence, X]fcLi k"^P{k) diverges 
for all m ^ 7 — 1. 
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2.5. Amenability, clustering, and small world properties 

The next property of our graphs which we are going to address is amenability. To introduce it 
we need one more notion. Let G = (V, E) be a countable graph with node set V and bond set E. 
For a finite A C V,hy OA we denote the set of nodes which are not in A but have neighbors in A. 
This set is the outer boundary of A, whereas the elements of A which are neighbors to the elements 
of dA constitute the inner boundary of Z\. As usual, by \A\ and \dA\ we denote the number of 
nodes in these sets. The graph G is said to be amenable if there exists a sequence of finite node 
sets {/4fe}fcgN , such that 

lim ^ = 0. (20) 

fe— > + CXD \Ak\ 

If such a limit is positive for any sequence {Ak}ken , the graph is called nonamenabile. The lattices 
Z'^, d ^ 1, are amenable graphs and for such sets one can take the cubes 

Ak^{i^{ii,...id)eZ'' : \ij\^Nk, j ^l,...,d}, 

such that iVfc-|.i > Nk and Nk — > +oo. In this case, \Ak\ ~ N^} and \dAk\ ~ -^^"^ and hence (1^ 
holds. Sometimes, sequences for which (|20p holds are called Van Howe sequences, see e.g. |21| . 
Cayley trees, except for Z, are nonamenable. Let us turn now to our graphs. Due to their hierarchical 
structure, it is convenient to check (|20p for the sequence of node sets of Afc , that is for {Vk}k€N ■ 
By the construction of Afc , the inner boundary of each Vk is exactly the set of all its external 
nodes, the number of which is equal to the number of nodes in the corresponding motif, i.e. it is 
q. By construction, one of them becomes an external node of Afc-|_i , and the remaining q — 1 ones 
become inner nodes of Afc_|_i . For all the motifs M, , j = 1, . . . , 5, we have the degree of any node 
therein being at most 3, see figure [T] Then for all our graphs, 

maxnfz) ^ 6k, 

c.f. (fT4|) . At the same time, \Vk\ ^ q^ , g = 3,4, see (|9]), which immediately yields that all our 
graphs are amenable. 

Now we study the clustering in our graphs. For non-random graphs, the clustering coefficient is 
defined as follows. For a given node i E V of degree n{i), let N{i) be the number of bonds linking its 
neighbors, which is the number of triangles with vertex i. Clearly, N{i) ^ n{i)[n{i) — l]/2 and the 
maximum value of this parameter is attained for complete graphs where each node is a neighbor 
to any other one. Thus, the quantity 

Qi^) ^4^T^ ' (21) 
n{i)[n{t) - 1] 

characterizes clustering at node i. Then we define the clustering of our graphs as 

Q= lim pi^E^W- (22) 

Note that for many graphs, e.g., for trees or bipartite graphs, one has Q{i) — for any node i, 
see also [2^, ^23!|. For random graphs, the degree n{i), as well as the parameter N{i), are random. 
Since in this case the calculation of Q is much more involved, we postpone it to the future. Here 
we only compare the values of Q obtained for the bare graphs with those for fully decorated ones. 

For the graph based on Mi and a node i € V^''\ I = 1, . . . ,k — l, we have: (a) n{i) = 4, N(i) = 3 
for I = 1, and N{i) = 2 for I ^ 2; (b) n{i) = Al, N{i) = 41. Here (a) and (b) correspond to a 
bare graph and to a fully decorated graph, respectively. These numbers follow directly from the 
construction of the graphs. The number of elements in each v/''' is given in (|12p . which allows one 
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to compute 

(a) Q = , lim I ^ + + -^1^1 = 0.4444 . . . , (23) 

1 3^/"' + 1 

(b) Q ^ 2 ■ 3~^/^arctan3~^/^ ~ — In 0.525897. (24) 

3 ' 3 ' ~~ 1 

For the bare graph based on , we have N{i) = for all nodes; hence, Q = in this case. For 
the fully decorated graph based on M3 , we observe that two nodes of ¥2^^ have neighbors only in 
V2^\ and the remaining four nodes have two neighbors in ¥2^^ and two — in V2'^\ Since for the 
nodes i G V2^\ we have N{i) = 0, like in the case of the bare graph, we should consider these two 
groups separately. Thus, we split each vl^\ Z = l,...,fc — 1, into V^'q and V^''l , where the first set 
consists of the nodes which have neighbors only in V^f g . The elements of V^^l have also neighbors 
in V^'' ■* with V > I. The number of nodes in these sets are: 

For all i e V^^'\ ? = 1, . . . , fc - 1, we have n{i) = 4/. At the same time, N{i) = 4(/ - 1) for i G vl:^l , 
and N(i) = 1 + 4(/ — 1) for i G vl^\ . Putting all these numbers together we obtain 

For the bare graph based on M5 one can obtain for i £ Vj!'\ ^ = 1, 2, . . . , fc — 1: n(i) = 6, Nii) — 8 
for I = 1, and n{i) — 6, N{i) = 6 for ^ ^ 2. For the fully decorated graph based on M5 we have for 
all internal nodes n{i) — 61, and N{i) = 9 for I = 1, and N{i) = 121 + 7 for I ^ 2. Hence, we obtain 

(b) Q « 0.554145. 

There exists one more property of real world networks, which Erdos-Renyi type graph does not 



share, see e.g. |20l. 124] . It is the so-called small-world property, illustrated by Pal Erdos himself 
in the following way. All authors of mathematical papers are given the Erdos index which for his 
coauthors is equal to one. For their coauthors who are not Erdos' coauthors, this index is equal 
to 2, and so on. Everyone whose papers are indexed by the Mathematical Reviews can check his 



own index at the web-site http:/ /www. ams.org/mathscinet/. It turns out that for many authors it 



ranges from 2 to 6 (for the second named author of this paper it is 4) . To formulate the small-world 
property one needs the following notion. A path in the graph is a sequence of nodes such that every 
two consecutive elements of this sequence are neighbors to each other. The length of the path is 
the number of such consecutive pairs, which is equal to the number of bonds one passes on the 
way from the origin to the terminus. If every two nodes can be connected by a path, the graph 
is said to be connected. For the given two nodes, i and j, the length of the shortest path which 
connects them is said to be the distance p{i,j) between these nodes. Informally speaking, a graph 
G — {V,E) has the small-world property (or is a small-world graph) if every two nodes i,j G V 
are "not too far" from each other. More precisely this property is formulated as follows. An infinite 
graph G has a small-world property if there exists a sequence of its connected finite subgraphs 
{GfcjfcgN with the following property. Let diam(Gfe) = maxij-gy^. p(i,j) be the diameter of Gk, 
fc e N, and (rik) be the average value of the node degree in Gk , that is, (uk) = 2\Ek\/\Vk\. Then 
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the sequence {GfcjfegN , and hence the graph G, are said to have the small-world property if there 
exists a positive constant C such that for all fc G N, 

diam(Gfe)<Clog^„^)|Vfc|. (27) 

This means that in such graphs, the distances between the nodes scale at most logarithmically 
with the size of the graph. For our graphs, the results on this item are given in the last three rows 
of table [U The diameters of presented therein were calculated in a routine way for the cases 
where the graphs are bare (p = 0) or are fully decorated {p = 1). In the former case, neither of our 
graphs has the small-world property. However, this property holds for all fully decorated graphs. 



Table 1 . The structural characteristics of the families of hierarchical graphs. 



motif 


Afi 


Ma 


Ah 


M4 


M5 


\Vk\ 


1(3*^-1 + 1) 


2(4*^-1 + 1) 


2(4*^-1 + 1) 


2(4'^-! + 1) 


2(4*^-1 + 1) 


\Ek\ 


1(3^-1 - l)p 


|(4'=-i - l)p 


1(4^-1 - l)p 


5.4^-1 
|(4fe-i _ i)p 


6. 4*^-1 
2(4'=-! - l)p 


(nk) 


4: + 2p 

A l+P 


4+Ip 

4 3+2p 


a + Ip 

4 3+2p 


5+Ip 

5 3+2p 


6 + 2p 

3+2p 




3 4'=-i + l 


3 4''-i + l 


3 4''-i + l 


^4fc-l+l 


diam(Afe), p = 


2k-i 


2k 




2k 


2k-i 


diam(Afe), p = I 


k 


k + 1 


2{k - 1) 


k 


k 


C,p = l 


(logeS)-! 


(log6 4)-i 


2(log6 4)-i 


(log7 4)-i 


(log8 4)-i 



In table [T] we collected the structural characteristics of our graphs based on the motifs Mi , . . . , 
M5 given in figure [T] Its first three rows are obtained from the formulas (O and (fTT)) . 

3. Graph structure and the Ising model properties 

As was mentioned above, there exists a profound connection between the properties of Gibbs 
random fields of models based on graphs and the structural properties of these gra phfl For the 
hierarchical lattices, the notion of the Gibbs random field of the Ising model was introduced 
in [l5l|. In the physical terminology, each (pure) Gibbs random field corresponds to a state of 
thermal equilibrium of the model, see j25| for more details. Accordingly, the existence of multiple 
Gibbs random fields corresponds to the existence of multiple equilibrium states and hence of phase 
transitions. If there is no interaction between the spins, the Gibbs random field is unique. However, 
if the interaction is strong enough and if it is effectively propagated by the underlying graph (high 
enough "connectivity"), the Gibbs fields can be multiple. The condition of high connectivity is 
essential, which can be seen from the fact that the Ising model on the one-dimensional lattice Z, 
considered as a graph with the natural adjacency relation, bears only one Gibbs field and hence 
no phase transitions can occur in this case - no matter how strong the interaction is. With regard 
to these arguments, we can divide all graphs into three groups according to the following property 
of the Ising model defined on these graphs: 

''Of course, we are speaking now about countably infinite graphs. 
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(a) There exists only one Gibbs random field for all interactions. 

(b) There exists only one Gibbs random field for weak interactions, and multiple Gibbs 
random fields for strong interactions. 

(c) There exist multiple Gibbs random fields for all nonzero interactions. 

As was mentioned above, the lattice Z belongs to the first group. The lattices with d ^ 2, as 
well as all Cayley tree^ belong to the second group. An example of the graph belonging to the 
third group can be found in section 4 of [26l . Clearly, the above classification should be made more 
precise for random graphs, since in this case the existence of a Gibbs random field is a random 
event. We rettirn to this issue below. 

The Ising model defined on the graphs from the second group may have a critical point, which 
separates the regime of uniqueness (weak interactions) from the regime of non-uniqueness (strong 
interactions). At their critical point, the Gibbs random fields have "unusual" properties, which 
can be detected without explicit construction of these fields. One of such properties is the so- 
called self- similarity, which manifests itself in the appearance of unstable fixed points of some 
(renormalization) transformations. It turns out that the hierarchical graphs are quite suitable for 
such a study — that it why they appear in this context. Thus, one can have the following criterion: 
if the Ising model has a critical point, the graph belongs to the second group. If not, then the 
graph belongs either to the first or to the third group. The latter cases can be distinguished by an 
additional study. 

In contrast to the hierarchical lattices introduced and studied in fl6j . our bare graphs belong 
to the first group. This can be seen from the analysis which we present below — a direct proof 
of such a statement will be done in our forthcoming publication. Thus, the interaction along the 
decorating bonds plays the main role in the possible appearance of phase transitions in the Ising 
model defined on our graphs. In view of this fact, we will assume that the solid bonds bear the 
interaction K, whereas the dotted bonds bear the interaction L. As the bonds of the latter kind 
appear at random, we can think of the corresponding model as of the one defined on the fully 
decorated graph with random interactions along the decorating bonds. Thus, in this model, the 
interaction along each bond {i,j) is a random variable, denoted by L^- , which takes values and 
some nonzero L gM. with probabilities l—p and p, respectively. For different bonds, these random 
variables are independent. 

For a given k G N, the Ising model on the fully decorated graph Ak is defined in the usual way 
by assigning spin variables cr^ = ±1 to the nodes i G Vk and by setting the Hamiltonian 

- l3Hk = hY^a,+K a,(jj + ^ iyCr.aj , /i, e R, fc e N, (28) 

where h is an external field. In accordance with the above arguments, the third summand cor- 
responds to the interaction along decorating bonds. We have included the first term in view of 
the following arguments. As is known, the Ising model on the lattices Z"^, d ^ 2, exhibits phase 
transitions only ii h — 0. This is also the case if the graph is amenable and quasi-transitive. For 
nonamenable graphs, this model may have a phase transition for nonzero h, see |27|. 

By the construction of our graphs, each has q external nodes, ii, . . . ,iq , and the correspond- 
ing number of internal nodes, see table [T] Such external nodes can be considered as a boundary of 
Afc . Correspondingly, we call the spin variables external or boundary (respectively, internal) spins if 
they are assigned to external (respectively, internal) nodes. For a fixed configuration of the external 
spins (Tij , . . . , (Ti^ , we consider 

Z^{a,,,...,a,J^ ^ expi-pUk). (29) 

internal spins of 

This is the conditional partition function of the spin system in Afe , conditioned by the fixed 
configuration of spins on the boundary of A^ . Of course, it depends on the graph, i.e. on the motif 

^Except for Z, see page 247 in \2^ . 
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used in its construction. According to the main principle of the theory of Gibbs random fields 
applied in our context, see [25], the fact that such a field is unique is ensured by the conditional 
partition function asymptotic independence, as fc — > +00, of the boundary spins. For our graphs, 
Z'^ are random and hence also depend on lo. In this situation, one can apply the following two 
approaches: the above mentioned independence holds (i) for almost all uj; (ii) in average. The latter 
corresponds to the fact that the expected values of Z'^ , which we denote by (Z'j^), are independent of 
the boundary spins. These two approaches are called quenched and annealed disorders, respectively, 
see [28|, page 99. We will work in the annealed approach. It can be shown that the summands in 
are positive and separated from zero. Since the number of these summands increases to +00 as 
k — > +00, we have {Z'^) — > +00. Therefore, a weaker form of the uniqueness condition can be 

lim (ZjfK,...,a,J)/(Z,"(a,:^,...,a: ))=c, (30) 
fc— ^+00 ^ 

which holds for any two configurations of the boundary spins and for some c > 0, which depends 
on these configurations. The convergence in ([50)) should be stable with respect to small changes of 
h and K, i.e. of the starting element Zi . In this case, the annealed limiting free energy per spin is 
independent of the boundary configuration. Correspondingly, a necessary condition for the latter 
quantity to be dependent on the boundary spins is that 

lim {Z^{a.,,...,a0)/{Z^K,---,<JI))^+^, (31) 

for some ct^j , . . . , and cr^^ , . . . , cr^^ . We say that the Ising model defined on our graph has a 
critical point, if there exists c' > 0, distinct from c in (|30p . such that for K = K^, and h — h^, , 

lim (Z,^(a,,,...,a.J)/(Zfe"«,...,a,:j)=c'. (32) 

Here if* and are certain values of the interaction intensity and the external field. The conver- 
gence in (|32p should be such that for arbitrarily small deviations from the point {K^, , /?,*) one has 
either (|30p or ([HT]) instead of ((5^ . Note that we do not exclude that /i* 7^ since our graphs are 
not quasi-transitive and hence the result of j23| may not be applicable in this case. 

4. Critical points of the Ising model on Mi based graph 

For k ~ 1, we have only basic bonds and no inner nodes. Thus, 

Zi{a, b, c) = exp [-;9'Hi(a, b, c)] = exp[i4:(a6 + ac + be) + h{a + b + c)], (33) 

where we use the same notations for the nodes a, b, c, see the left-hand graph in figure |5] as well 
as for the corresponding spin variable, that is, for short we write a = aa , b = ab , ect. For k = 2, 
we have, see the middle graph in figure |2] 

— /?'H2(a, b, c, a, fi, 7) = — /?'Hi(a, 7, (3) — /3Hi(7, ^7 ~ P'^iiP, ot, c) -f L'^b^^ + L^Jic + L^^c^-c. 
Thus, to obtain ^2(0, 5, c) we have to sum out the internal spins, see ((29)) . That is, 

Z^ (a, 6, c) = exp [L^^^ab -f L^^ac + L^hc] ^ Zi (a, 7, /5)Zi (7, ^ a)Z^{l3,a,c), 

where the sum is taken over a, /3, 7 = ±1. According to the hierarchical structure of the underlying 
graph, for all fc ^ 2 the graph A^ is obtained from the subgraphs Afe_i exactly in the same way. 
This yields the following recurrence for the partition functions ()29p 

Zk{a,b,c) = R%{a,b,c) ^ Zfc_i(a, 7, /3)Zfe_i (7, 6, a)Zfe_i(^, a, c), (34) 
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where 

(a, &, c) = exp [L^,a6 + L^.ac + 5c] , (35) 

and ^i(a, 6, c) is given in ([55]) . From the independence of for different bonds, it follows that the 
random variables Z'^_^ on the right-hand side of p4p are independent, which allows us to compute 
the expectations and obtain 

{ZtiaAc)) - {RtiaAc)) {Zt^MP.l)){Zt-i{^Al)){Zt-i{oc,p,c)). (36) 

We recall that the random variables L^- take the values L 7^ and with probabilities p and 
1— p, respectively. In view of the model symmetry, c.f. (|33p . the system p6p involves the following 
variables 

Ak = (^,"(1,1,1)), B,^{Z^{l,l,-l)), (37) 
Ck - (Z,"(-l,-l,l)), Z?fc = (Z-(-l,~l,-l)), 
which have the initial values 

Ai-eWM^ i?i=e-^+\ Ci-e-^'-'\ i^i = e''*'^-''). (38) 
By direct calculations, 

(i?^ (1,1,1)) = ipcxp{L) + l-p)\ (39) 
(i?^(-l,l,l)) = (pexp(L) + l-p)(pexp(-L) + l-p)2. 
Thereafter, the recursion p6p can be written in the form 

a^fc+i = ^7^7 7> (40) 



2/fe+i 
2fe+i 

with the initial conditions 



Py{Xk,Vk,Zk) 

Q{xk,yk,zk) 

^ Pz{yk,zk) 
Q{xk,yk,Zk) ' 



XI -eW), y^^,2H^ z,=e^^-'\ (41) 
where we have used the notations 

Xk=Ak/Ck, Vk^Bk/Ck, Zk = Dk/Ck, fceN, (42) 

and 

.= fi4±i^V. (43) 
The polynomials in (|40p have the following form 

PAx,y) = + 3x1/2 + 31/2 + 1, (44) 
Py{x,y,z) — x'^y + + 2xy + y'^z + 2y + z, 

Pz{y,z) = y^ + 3y + 3z + z^, 
Q{x,y,z) = xy'^ + x + 2y^ + 2yz + z'^ + 1. 
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For every t > 0, the mapping (xfc , jjk , Zk) ^ (xk+i , yk+i , Zfe+i) = Tt{xk , yk , Zk) , defined in (gO]) 
and (|44p . maps the first octant of into itself. The starting point {xi ,yi ,zi) lies on the surface 
S defined by the equation, c.f. (|¥T|) . 

X = y^z, (45) 

which, however, is not preserved by this mapping, which can be checked directly. Let 7t be the 
invariant set of the map Tt , that is, the set of points {x,y,z) such that Tt{x,y,z) — {x,y,z). 
Suppose that the intersection of Tt with S is non-void. If (x* , , z^) belongs to this intersection 
and {xi , yi , zi) = (a;* , , 2*) , then {xk , yk , Zk) — (x* , , z*) for all fc G N. Let us find such 
points. From the second equation in (|40p for x = y^z we get 

yQ{y^z, y, z) = Py{y'^z, y, z), 

which can be transformed into the following 

+ y^z){l ~ yz) = {y + z){l- yz). (46) 

A solution of the latter equation is z = thus x = y"^ , which being used in the first equation 
in (pOj) yields (1 + y^)^ = t{l + y^)^ and hence t — 1. For p > 0, the latter imphes L = 0, and also 
K ^ 0, see (|4T|) . This fixed point corresponds to a noninteracting system of spins in an external 
field and thereby has nothing to do with critical points. For z ^ 1/y, from we get 

zil-y^)iy^ + y' + l) = -yil^y^), (47) 

which for positive y and z has a unique solution y = I, and hence x — z. We insert this into the 
first equation in (|40p and obtain 

x^ + 3x + 4 x'^-x + A 

For i = 1, the only solution is x — 1, which corresponds to L — K = {). The choice t — 1 also 
corresponds to a bare graph (p = 0); hence, the bare graph based on Mi has no critical points 
and no phase transitions. As we have already mentioned above, in a separate work we prove that 
the Gibbs states of the Ising model based on this graph are always unique. For t <1, that is, for 
L < 0, (P5]) has only one positive solution 



-3 - t + \/9 + 22* - 15*2 
X. = < 1, (49) 

which corresponds to if* < 0. It is stable and {xk ,yk , Zk) ~> (x*,l,a;*) as — > +00, for any 
{xi ,2/1 ,zi) S S. Thus, for L < Q, the Ising model has no critical points and no phase transitions, 
even for the fully decorated graph. This can be caused by a frustration, which takes place in an 
antiferromagnetic Ising model on triangles. For t > 1, figure [S] presents the graphical solutions 
of (|48p . They have the following form 



(1) _ 3 + * - V9 + 22* - 15*2 _ 3 + * + V9 + 22* - 15*^ 

~ 2(*^T) ' " 2(*^T) • 

By direct calculations, we get (l)[{xi^^) < 1 and i^ti^^^^) > Ij see figure |5] Hence, xi^^ is stable, 
whereas xl^-* is unstable. Note also that xl^'' — t- 1 and x^^^ — )■ 00 as * — 1. The solutions ((50)) exist 
and are considered to be distinct, provided * G (1,9/5). This yields the condition 

^ ^ -\/5exp(L) — 3exp(— L) + 3 — a/5 ^ ^ 

Since ^(L) is a decreasing function, the equation 'ijj{L) = 1 has only one solution 

= iln^. (52) 
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Then, for L ^ L» , the critical point exists for all p E (0, 1]. For such L and p, let K^,{L,p) be the 
solution of the equation 

eM'iK)=x'-/K (53) 
if xi < exp[4if,(i,p)], 

if xi^e^p[AK^{L,p)], (54) 
if xi > exp[4ii',(L,p)]. 

The conditions in (|54p can be formulated directly for K, e.g. Xk — xi^^ if K < K^,{L,p). Thus, 
zl^'' is the so-called high-temperature fixed point. For L > L^t , we have ip{L) < 1. Hence, the 
condition (|5T|) is not satisfied for p G {iJ^{L), 1], which also includes the fully decorated graph with 
p — 1. For such L and p, the whole graph of 0t lies above the line (l)t{x) — x, which means that 
Xk —7- +00 for all initial xi ^ 0. Thus, for L > and p G (ip^L^), 1], the model has no critical 
points and the Gibbs random fields are multiple for all K G M.. For p — ip{L^:), we have t — 9/5 
and a;!^"* = a;!^' = 3, which corresponds to A', — (ln3)/4. In this case, Xfc — >■ 3 if AT ^ if, , and 
Xk — > +0O ii K > . 



Then 



X 



(1) 



lim Xk = 



5. Concluding remarks 

5.1. Construction and structural properties 

The graphs introduced in this paper are constructed from motifs Mi , . . . , in an algorithmic 
way, like the hierarchical lattices introduced in fisj , and employed in |Jj£] where they were supplied 
with long-range (decorating) bonds. However, our main construction principle — to replace nodes 
with the graphs of the previous level — differs from the one used in those papers, where the graphs 
of the previous level replaced the bonds. Hierarchical graphs of this kind model the real world 
networks with the so-called modular structure, which are organized in tightly knit communities 
with relatively sparse connections between them. For more details herein, we refer the reader 
to [i^^O] and to the publications cited therein. In contrast to the graphs introduced in [l6|, see 
also |29|, our decorated graphs are not scale-free. For them, the node degree distribution is of 
intermediate type such that all the moments exist. The corresponding characteristic function can 
be extended to a meromorphic function, analytic in the complex neighborhood of zero, see psp . 
whereas for scale-free graphs such functions are nonanalytic at zero. All our graphs are amenable, 
which correlates with the property just discussed, see (fT4)) . ([9]), and (fl9|) . An interesting property 
is that, like for the graphs studied in [l^, adding decorations forms the graphs possessing the 
small-world property. 
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5.2. Critical points 

The appearance of critical points, and hence of phase transitions, for the Ising model confirms 
good communicating properties of the underlying graphs. For our graphs, in contrast to the hierar- 
chical lattices for which phase transitions take place even without decorations, critical points can 
exist only if the decorations are applied, with any p > 0. The most interesting related fact is that 
for the graph based on Mi , the antiferromagnetic Ising model has no phase transitions for any L 
and p. This is due to the generic frustration of such models — the elementary factor Rk{a,b,c) 
which appears in psp with negative L cannot be maximized. We expect that the same is true for 
the graphs based on . To check if this conjecture of ours is true we plan to study a Potts model 
on the same graph and q ^ 3, for which such a frustration is no longer actual. We also plan to 
look for critical points with nonzero h, as well as for critical points of the Ising model with the 
underlying graphs based on the remaining motifs. However, the corresponding problems are much 
tougher, so that the appropriate numerical means should be applied. 
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lepapxiHHi Bi/ina^KOBi rpac^n 3 MOTi/ieiB: CTpyKTypni B/iacTi/iBOcri ra 

Kpi/ITMHHi TOHKM MOfle/li IsiHTa 

M. KoTopoBi/iH, hOpiw K03l/mbKI/lPi 

iHCTMTyX MaXeMaTMKM, VHiBepCMTeT MapiT Kropi-CK/lOflOBCbKoT, f\\o6f\\H 

BBOflMTbca i BMBHaeTbca KJiac BunaflKOBi/ix rpacjjiB, 36yflOBaHi/ix b ajiropnTMiHHwti cnoci6 3 n'axw moti/ibIb, 
BHaCifleHMX y [Milo R., Shen-Orr S., Itzkovitz S., Kashtan N., Chklovskii D., Alon U., Science, 2002, 298, 
824]. KoHCTpyKL4iMHa cxewia Haraflye cxeMy, sacTocoBaHy y [Hinczewski M., A. Nihat Berker, Phys. Rev. E, 
2005, 73, 056125], sriflHO 3 aKoio KopoTKOC5i>KHi pe5pa e HeBi/inaflKOBi, xofli aK flOBroc5i>KHi pe5pa bwhw- 
KaKDTb He3ajie>KHO i 3 oflHaKOBOKD tiMOBipHicTKD. Oni/icaHO pjifl CTpyKxypHMX BJiacTMBOCTew rpa4)iB, cepefl 
!RKMX e po3nofli/i CTyneHiB, K/iacxepHicTb, aMeHaGi/ibHicxb, B/iacxuBicTb ricHoro CBiry. fl/in oflHoro 3 MOTi/iBiB 
BMBHacTbCiR KpMTMHHa TOHKa Moflejii bihra, BH3HaHeHoT Ha BiflnoBiflHOMy rpactii. 

K/iiOHOBi c/iOBa: ameHaGinbHicTb, posnofliji cryneHiB, KJiacrepHicTb, rpacfi ricHoro CBiry, MOftent hiHra, 

KpHTUHHa TOHKa 
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